:arXiv:hep-th/9605068vl 9 May"l996i 



THE PATH INTEGRAL FOR THE LOOP REPRESENTATION OF LATTICE 

GAUGE THEORIES 

J. M. Aroca 

Departament de Matematiques, Universitat Politecnica de Catalunya , Gran Capita, s/n Mod C-3 Campus Nord, 08034 

Barcelona, Spain. 

H. Fort and R. Gambini 
Institute de Fisica, Facultad de Ciencias, Tristan Narvaja 1674 ,11200 Montevideo, Uruguay 

(February 1, 2008) 

We show how the Hamiltonian lattice Zoop representation can be cast straightforwardly in the 
path integral formalism. The procedure is general for any gauge theory. Here we present in detail the 
simplest case: pure compact QED. The lattice loop path integral approach allows to knit together 
the power of statistical algorithms with the transparency of the gauge invariant loop description. 
The results produced by numerical simulations with the loop classical action for different lattice 
models are discused. We also analyze the lattice path integral in terms of loops for the non-Abelian 
theory. 
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I. INTRODUCTION 



The loop approach to abehan quantum gauge theories was introduced in the early eighties Later it was 
generalized to the non-abelian Yang-Mills gauge theory 1^. This Hamiltonian method allows to formulate gauge 
theories in terms of their natural physical excitations: the loops. The original aim of this general description of gauge 
theories was to avoid gauge redundancy working directly in the space of the gauge invariant excitations. However, 
soon it was realized that the loop formalism goes far beyond of a simple gauge invariant description. The introduction 
by Ashtekar |^ of a new set of variables that cast general relativity in the same language as gauge theories allowed 
to apply loop techniques as a natural non-perturbative description of Einstein's theory. In particular, the loop 
representation appeared as the most appealing application of the loop techniques to this problem [Q . The covariant 
worldsheet formulation of quantum gravity corresponding to the canonical loop representation is still unknow. This 
version of loop quantization would make the general 4-dimcnsional diffcomorphism invariance manifest and it will be 
probable more suitable to tackle some crucial problems as the black hole dynamics. Some progress in that direction 
has been reported recently 

The Hamiltonian analytical computation techniques for gauge theories have been developed during the last decade 
and they provide qualitatively good results for several lattice models On the other hand a Lagrangian approach 

in terms of loops has been elusive, due mainly to the non-canonical character of the loop algebra. This feature forbids 
the possibility of performing a Legendre transformation as a straightforward way to obtain the Lagrangian from 
the Hamiltonian. A path integral loop formulation will allow to employ the more powerful statistical computation 
techniques. Gauge invariant actions corresponding to the Villain form of the U(l) model were proposed recently |^ 
and generalized to include matter fields | |ic| ]. 

These actions has been used as a computational tool. A Metropolis Monte-Carlo algorithm was implemented fixing 
the acceptance ratio, as it is usual in random surfaces analysis. The U(l) model was studied for different lattice sizes 
1^, imposing periodic boundary conditions. Simple thermal cycles showed the presence of a phase transition in the 
neighborhood of Pv = 0.639. In the case of of U(l)-Higgs model, the (3-^ phase diagram was mapped out jlOj. In 
both models, by virtue of the gauge invariance of this description, the equilibrium configurations were reached faster 
than with the ordinary gauge-variant descriptions. Furthermore, the absence of the typical strong metastability of 
Monte Carlo analysis , which makes very difficult the numerical analysis of the phase transition critical exponents, 
was observed. This fact is one of the main advantages for using the loop action. 

In this paper we want to address the issue of the explicit correspondence between these actions and the Hamiltonian 
loop representation using the transfer matrix techniques. 

This paper is organized as follows. In section 2 we show how the loop description, originally introduced in the 
Hamiltonian formalism, can be cast in the lattice path integral formalism. We illustrate this in detail for the Abelian 
case (compact electromagnetism) . The path integral of lattice U{1) theory is expressed as a sum of the world sheets 
of electric loops. We discuss the connection of this classical loop action with the Nambu string action. In section 3 
we consider the extension to the path integral loop formalism to the case of non- Abelian Yang-Mills fields. In section 
4 we conclude with some remarks. 



II. THE LATTICE PATH INTEGRAL IN TERMS OF LOOPS: ELECTROMAGNETISM. 

The loop based approach of ref. describes the quantum electrodynamics in terms of the gauge invariant holonomy 
(Wilson loop) 

M/(7) = exp[ze/A,(y)dy°], (1) 

where 7 is a spatial loop at constant time t. E°'{x) is the conjugate electric field operator. They obey the commutation 
relations 



[i?"(x), iy(7)] 5{x- y)dy^'W{-i). (2) 

These operators act on a state space of abelian loops V'(7) that may be expressed in terms of the transform 
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1/^(7) = / d^[A] < 7 I A X A I / d^[A]',p[A] exp[-ze i Aady'']. (3) 



This loop representation has many appeahng features: i) It aUows to do away with the first class constraints of 
gauge theories i.e. the Gauss law is automatically satisfied, ii) The formalism only involves gauge invariant objects 
i.e. no gauge redundancy, iii) All the gauge invariant operators have a transparent geometrical meaning when they 
are realized in the loop space. 

When this loop representation is implemented in the lattice equation (j^) becomes 

[4,M/(7)] = 7V,(7)#(7), (4) 

where i denotes the links of the lattice, Eg the lattice electric field operator, W{'^) — W^^^ U{£) and Ni{'j) is the 
number of times that the link i appears in the closed path 7. 

In this loop representation, the Wilson loop acts as the loop creation operator: 

W^(7') I 7 >=| 7' .7 > . (5) 
The physical meaning of an abelian loop may be deduced from ^ and (^, in fact 

Ei I 7 Ni{-f) I 7 >, (6) 

which implies that | 7 > is an eigenstate of the electric field. The corresponding eigenvalue is different from zero if 
the link I belongs to 7. Thus 7 represents a confined line of electric fiux. 

The U{1) lattice Hamiltonian can be written in terms of both previous fundamental operators as 



Now, starting from the path integral for the Wilson U(l) lattice action, let us show how to set up the loop path 
integral. 



Zw = J [dee]exp{^J2 



cos 6'p), (8) 



where the subscripts £ and p stand for the lattice links and plaquettes respectively, (3 = (e is the electric charge of 
the electron) and 9i is a compact variable G [— tt, tt] attached to links and Op is its "discrete curl". Fourier expanding 
the exp[/?/2 cos6'] we get 



p rip 

where the /„ are modified Bessel functions. Zw can be written as 



^[de.]^E^".(/^)e"''^ (9) 

•J — TV ^ 



/TT 
[d9e]exp{Y^ ln/„,(/3))e' 

inpi P 



<n,'V9t> 



(10) 



where we use the notations of the calculus of differential forms on the lattice. In the above expression: 6* is a 
real compact 1-form defined in each link of the lattice V is the co-boundary operator -which maps fc-forms into 
(k + l)-forms-, np are integer 2-forms defined at the lattice plaquettes and consider the scalar product of p- forms 
< / I 5 >— X^Cfc /(c)5(c) where the sum runs over the /c-cells of the lattice (cq sites, ci links and so on ). Under 
this product the V operator is adjoint to the border operator S which maps k-forms onto (k-l)-forms and which 
corresponds to minus times the usual divergence operator i.e. 

<f\Vg>^<Sf\g>, (11) 
< V/ I 5 >-< / I <5g > . (12) 
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Using equation (|12|) and integrating over 9^ we get a 5{5np). Then, we arrive to the following functional of the integer 
2-fornis n 

Zw[n]^ exp(^ln/„^(/3)), (13) 

{np;i5nj,=0} P 

the constraint Sup — means that the sum is restricted to closed 2-forms. Thus, the sum runs over collections of 
plaquettes constituting closed surfaces. 

An alternative and more easy to handle lattice action than the Wilson form is the Villain form. The partition 
function of that form is given by 

Zv = y'M0]^exp(-^ II V0-27rs|n, (14) 
where || ... |P=< ....... >. If we use the Poisson summation formula 

/oo 
d0/(0)e2"^" 

and we integrate the continuum variables we get 

Zv = (27r/3y)-^''/2 / [dO] y exp(--^ <n,n> +i <n,Ve >), (15) 
J V 

where Np in the number of plaquettes of the lattice. Again, we can use the equality: < n, V0 >=< Sn^O > and 
integrating over 9 we obtain a 5{5n). Then, we get 

Zv[n]^{2n(3v)-''-'^ exp(-^ <n,n>). (16) 

Let us first discuss the (2+1) dimensional case. In figure |^ we show the result of intersecting one of the previous 
surfaces with & t — constant plane. We can get: i) a spatial closed path or loop (Figure l-(a) ) or ii) an open surface 
connecting the loop 7^, obtained by translating the loop 7t_i by one temporal lattice unit, and 74 (the shaded surface 
in Figure l-(b) ). A loop 7* living on the t slice is properly specified by the temporal plaquettes which leave this slice. 
This is equivalent to say that in two spatial dimensions, given the loops at time t and time i + oq, the surface which 
connect them is unambiguously defined. 
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In more than 3 space-time dimensions the situation is different, the loops ^t-i and 74 do not define a unique world 
sheet connecting them. Thus, if we consider the intersection of one of the world-surfaces with a t = constant plane 
we can get: i) a loop 74 ii) an open surface connecting the loop 7^ and 74 iii) a closed spatial surface. The situation 
is completely analogous to the path integral of a free particle on a (2-1-1) dimensional lattice i.e. the world lines 
connecting points at different times intersected with a, t — constant plane give: i) points or ii) open paths connecting 
the intersecting point of the timelike link at pt-i with surface t — constant and the point pt and iii) loops associated 
with different choices of the previous mentioned open paths. i-i 

Now we will show that the expression ( p^ ) is a possible "loop path integral" Zl First, it is easy to prove that 
the creation operator of these loops is just the creation operator of the loop representation, namely the Wilson loop 
operator. Repeating the steps from equation (nj) to equation (Rq) we get for < W{-/t) > 



< W^(74) i(2^/3^)-^./2 J2 exp(-^<n,n>). (17) 

n 

{6n = 7t) 

This is a sum over all world sheets spanned on the loop jt- Secondly, by means of the transfer matrix method let 
us show that we re-obtain the Hamiltonian (|^) from (p^. As we wish to consider the continuous time limit of the 
previous lattice Euclidean space-time theory, we introduce a different lattice spacing ao for the time direction. The 
couplings on timelike and spacelike plaquettes are no longer equal in the action i.e. we have two coupling constants: 
/3o and (3s. The temporal coupling constant f3a decreases with whilst the spatial coupling constant Ps increases 
with flQ. We wish to find an operator T over the Hilbert space of loops {I7 >} such that, the loop path integral Zl 
is given by 

= n < '^t+ao\T\jt > ■ 
t 

T is related, when oq is small, with the Hamiltonian H by 

f (xe""°^+°(''«^. (18) 
The transfer matrix between times t and t + in the loop representation of kets I7 > can be written as 

<7t+ao|T|7t >=exp[^^7i2--i- ^ nl], (19) 

where St+ao and St are the surfaces enclosed by the loops ^t+ao and 74 respectively. Due to the fact that the surfaces 
are closed the integers Up of the temporal plaquettes which depart from the loop 74 at time t and arrive to the loop 
jt+ao at time t + ao are equal to the number of times the spatial link £t appears in the loop 74, Ng. Therefore, the 
equation (^9|) can be re-written as 

<7t+ao|r|7t>=exp[-^^iV2(74)-^ nl]. (20) 

The kets |74+ao > and 74 > are connected by 

l7t+ao>= n ^;'"l7t>- (21) 

p£AoS~i 



Using equations (||), (||) and (^ ( p5[ ) we get 



E n K^^M-^Xe!-±- nl], (22) 



* Another possible path integral in terms of loops is given by (p^. 
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where {AoS'^,} are all the possible modifications between the surfaces enclosed by the loops jt+ao and jt i-e. it denotes 
the set of possible configurations {rip} of integers attached to the plaquettes p G AqSj . Then, 

^ = nE[^;''^-P(-2F E nl)]eM-^T.^"]- (23) 



To obtain a proper classical limit we should take 



/3o = ^ (24) 
5 flo 



" 2 ln(25 Vao) ' ^^^^ 
where a continues to denote the spacelike spacing. This implies that for ao small the operator T is given by 

f = exp{-ao[ E + ^ E(^f + ] + 0(«o)}, (26) 

i.e. we recover the Hamiltonian (^. This definitely confirms that the expression of the partition function of compact 
electrodynamics in terms of the world sheets of loops: the loop (Lagrangian) representation. 

From equation (|6|) we can observe that the loop action is proportional to the quadratic area A2: 

5L = --^A2---5^^n2 = --^<n,n>, (27) 
Pv PV ^ Pv 

i.e. the sum of the squares of the multiplicities Up of plaquettes which constitute the loop's world sheet S. It is 
interesting to note the similarity of this action with the continuous Nambu acticm or its lattice version, the Weingarten 
action which are proportional to the area swept out by the bosonic string CJ. 

III. THE NON-ABELIAN LOOP ACTION 

Let us see how can be extended the path integral loop description for the non-Abelian case of Yang-Mills theory. 
The path integral for the Wilson action for a general non-Abelian compact gauge group G is given by 

Zw-- j [dUt] exp[/3 Re(tr Up) ] , (28) 



where the Ue € G and Up = Yiiep^^- Equation ( p8| ) reduces to (|^) for the case G = U(l). The analogous of 
the Fourier expansion for the non-Abelian case is the character expansion. The characters XriU) of the irreducible 
(unitary) representation r of dimension dr, defined as the traces of these representations, are an orthonormal basis 
for the class functions of the group i.e. [Q3[ 



dUxr{U)x:{U) ^ Sr. (29) 
Y,drXr{UV~')^S{U,V). (30) 

r 

In particular, as a useful consequence we have 

dr f dUXs{U)Xr{UV-') = SrsXAV). (31) 



^The relation between the surfaces of the Wilson action and those of Weingarten action has been analyzed by Kazakov et al 
in ref. IT2I. 
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By means of the character expansion we can express 

exp{/3^ Re[x(;7p) ]} = [] E ^rxAUp), (32) 



p p r 

with 



d[/x:(C/)exp(/3x(C/)). (33) 
For instance, in the case of G = SU(2) the gauge fields can be parametrized as 

U = cos ^6 + iaaUa sin ^6, < 9 < An , 

and the corresponding irreducible representations are classified by a non-negative integer or half-integer spin j i.e. 
r = j and the characters are given by 

„iU) ^ ,34, 
sm -^0 

A direct application of (^3|) yields the cj in terms of modified Bessel functions, and therefore we can express 

j2j,+iiP) sin(j-f i)0p 



Zjw — 



p jp 



= E Uim. + 1)^^^] / m n (35) 
fe} P p J p ^^^^ 

A given subset of plaquettes is homeomorphic to a simple surface if any link bounds at most two plaquettes of this 
subset. The links bounding exactly one plaquette make up the boundary of this surface. Any configuration can be 
decomposed as a set of maximal simple surfaces by cutting it along the links bounding more than two plaquettes. 
In principle, there are two possibilities for the boundary curves: either a free boundary, bounding only one simple 
surface or a singular branch line along which more than two simple surfaces meet. In fact, relation (|2^) forbids the 
existence of free boundaries for non trivial configurations contributing to the path integral. 

The integration over the internal links of the simple surfaces is performed using (pT|). Note that the plaquettes of 
a simple surface component should carry the same group representation. After integrating over all the inner links of 
the simple components one gets an expression involving only the links of the boundary i.e. something proportional to 



n (^boundary)- 
boundaries 

What follows is the integration of gauge fields along the singular branches which gives rise to the Clebsch-Gordan 
coefHcients coupling the different representations of the considered gauge group G. For instance, imagine that there is 
only one singular closed branch line which is the common boundary shared by n simple surfaces with representations 
ri , r2 , . . . , r„ . The integration over gauge fields produces a factor A'^ , . . . .r„ , which counts the number of times the trivial 
representation is contained in the product ri (g) r2 03 ... (8) r^. 

Cases in which different singular branch lines meet in a point are considered in ref. Q . Each point of intersection 
involves a Racah-Wigner symbol. Thus, one can see that the Hamiltonian formulation associated to this action will 
be given in terms of a spin network of colored loops fl^ . A rigorous proof of this fact using the transfer matrix 
technique is still required. This is not straightforward because one still need to develop the lattice Hamiltonian 
formulation of Yang-Mills theory in terms of spin networks. 

We can also generalize the Villain form of the action for any gauge group using the heat kernel action , : 

exp{(3SHK) = nE^'-^-(^p)e^P[^^'-(2)/^/3] (36) 

p r 

where Cr{2) is the quadratic Casimir invariant for the representation r. For G = SU(2) the heat kernel action reads 

sinf i + — ]6 

eMf3SHK)=ll J2 + % xpH(j + l)/2/3], (37) 

sin 

P j=0,l/2,.. 2 

while only integers values of j are used for the SO (3) group. 
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IV. CONCLUSIONS 



As it was mentioned, the loop space provides a common scenario for a non-local description of gauge theories and 
quantum gravity. The loop approach is no more exclusively Hamiltonian , its Lagrangian counterpart is now available. 
A path integral action for the Yang-Mills theory in terms of loop variables is very valuable because it combines the 
geometrical transparency and economy of the loop description with the versatility to perform calculus. We have 
presented the state of the art in that program, which still is an open issue. 

The path integral approach to quantum gravity has very appealing features. In particular it may provide a more 
suitable framework for the development of useful approximation schemes for the study of black hole physics and it 
may allow to analyze issues such as the computation of the probabilities for a change of the spatial topology that 
seem to be very difficult to formulate in the canonical approach. Even though the the connection of the canonical 
loop representation of quantum gravity and the path integral approach is still an open problem the determination of 
the explicit form of the loop actions in gauge theories is an important step in this direction. An important remark is 
that the lattice framework seem to be unavoidable in order to have well defined loop actions. In fact in spite of the 
similarities with the Nambu actions, the loop actions for gauge theories involve cuadratic surface elements that are 
not well defined in the continuum. 

Finally, concerning the lattice loop action as a computational tool, we already mentioned that the results produced 
by numerical simulations for different models are very encouragging. In ref. |^ , the loop action corresponding to 
Villain form of U{1) model was considered. The extension of the Lagrangian loop description in such a way to include 
matter fields was also simulated The lattice path integral of [/(l)-Higgs model is expressed as a sum over closed 
as much as open surfaces. These surfaces correspond to world sheets of loop- like pure electric flux excitations and open 
electric flux tubes carrying matter fields at their ends. This representation is connected y a duality transformation 
with the topological representation of the path integral (in terms of world sheets of Nielsen- Olesen strings [|8| both 
closed and open connecting pairs of magnetic monopoles). Simulating numerically the loop action corresponding 
to the Villain form, the two-coupling phase diagram of this model was mapped out. The gauge invariance of the 
loop description bears the advantages of economy in computational time and the absence of the strong metastability 
previously observed in the ordinary Monte Carlo analysis. 
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